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Geometric  design  of  an  integrated  thermoelectric  generation-cooling  system  is  performed  numerically 
using  a  finite  element  method.  In  the  system,  a  thermoelectric  cooler  (TEC)  is  powered  directly  by  a  ther¬ 
moelectric  generator  (TEG).  Two  different  boundary  conditions  in  association  with  the  effects  of  contact 
resistance  and  heat  convection  on  system  performance  are  taken  into  account.  The  results  suggest  that 
the  characteristics  of  system  performance  under  varying  TEG  length  are  significantly  different  from  those 
under  altering  TEC  length.  When  the  TEG  length  is  changed,  the  entire  behavior  of  system  performance 
depends  highly  on  the  boundary  conditions.  On  the  other  hand,  the  maximum  distributions  of  cooling 
power  and  coefficient  of  performance  (COP)  are  exhibited  when  the  TEC  length  is  altered,  whether  the 
hot  surface  of  TEG  is  given  by  a  fixed  temperature  or  heat  transfer  rate.  The  system  performance  will 
be  reduced  once  the  contact  resistance  and  heat  convection  are  considered.  When  the  lengths  of  TEG 
and  TEC  vary,  the  maximum  reduction  percentages  of  system  performance  are  12.45%  and  18.67%, 
respectively.  The  numerical  predictions  have  provided  a  useful  insight  into  the  design  of  integrated 
TEG-TEC  systems. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

In  recent  years,  the  problems  of  global  warming  and  air  pollu¬ 
tion  have  greatly  stimulated  the  development  of  green  energy 
technologies.  In  these  technologies,  thermoelectric  devices  are 
receiving  a  great  deal  of  attention  duo  to  their  numerous  attractive 
merits.  For  example,  they  are  compact  and  environmentally 
friendly,  and  can  be  operated  easily  with  long  life  and  low  mainte¬ 
nance  cost  [1,2],  According  to  thermoelectric  effects,  thermoelec¬ 
tric  devices  have  two  different  operating  modes.  One  is  the 
generation  mode  using  thermoelectric  generators  (TEGs)  to 
directly  convert  thermal  energy  into  electrical  energy  [3],  The 
energy  source  of  TEGs  can  come  from  waste  heat  which  is  exten¬ 
sively  available  in  industry  [4,5],  The  other  is  the  cooling  mode 
whereby  thermoelectric  coolers  (TECs)  can  cause  a  temperature 
difference  for  cooling  applications  by  inputting  electrical  energy 
[6],  The  development  of  TEG  and  TEC  has  been  considered  as  a 
low-carbon  and  green  energy  technology,  and  they  have  been  suc¬ 
cessfully  applied  in  military,  aerospace,  and  industry  [7], 

In  general,  TEGs  and  TECs  are  individually  utilized  for  practical 
applications,  such  as  waste  heat  recovery  [8,9]  and  ceiling  cooling 
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system  [10],  Recently,  an  integrated  TEG-TEC  system  has  been  pro¬ 
posed  by  Chen  et  al.  [  1 1  ].  In  their  integrated  system,  TECs  were  dri¬ 
ven  by  TEGs  so  that  no  additional  electrical  power  source  was 
required  for  the  TECs.  They  provided  the  optimal  number  ratios 
between  the  TEGs  and  the  TECs  under  various  operating  conditions 
through  an  analytical  analysis.  Later,  several  studies  also  investi¬ 
gated  the  integrated  systems.  For  example,  Khattab  and  Shenawy 
[12]  constructed  an  experimental  system  and  successfully  used 
solar  TEGs  to  drive  a  TEC  all  year  round.  Meng  et  al.  [13]  used  an 
analytical  method  to  analyze  the  performance  of  a  thermoelectric 
heat  pump  driven  by  a  TEG.  They  found  that  the  heat  source  tem¬ 
perature  of  the  TEG  had  a  greater  effect  on  heating  load  than  on  the 
coefficient  of  performance  (COP).  On  the  contrary,  the  heat  sink 
temperature  of  the  thermoelectric  heat  pump  had  a  more  signifi¬ 
cant  effect  on  the  COP  than  on  the  heating  load.  Meng  et  al.  [14] 
also  adopted  another  analytical  method  to  investigate  the  influ¬ 
ence  of  physical  dimensions  of  thermoelectric  elements  on  the  per¬ 
formance  of  an  integrated  TEG-TEC  system.  Their  results  indicated 
that  performance  improvements  could  be  achieved  by  optimizing 
the  physical  dimensions  of  thermoelectric  elements. 

Geometric  design  is  a  crucial  issue  in  optimizing  performance  of 
a  thermoelectric  system.  For  a  given  thermoelectric  module,  the 
geometrical  optimization  of  the  thermoelectric  element  is  a  feasi¬ 
ble  way  to  maximize  its  performance  or  efficiency  [15],  Regarding 
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Nomenclature 

A 

cross-sectional  area  of  the  collector  (mm2) 

Te 

vector  of  nodal  temperatures  (°C) 

Ac 

foo 

environment  temperature  (°C) 

Cg 

COP 

concentration  ratio  of  solar  thermoelectric  generator 
coefficient  of  performance  of  the  integrated  system 

W 

width  of  thermoelectric  element  (mm) 

a 

depth  of  thermoelectric  element  (mm) 

Greek 

letters 

E 

electric  field  intensity  vector  (V m  ') 

a 

Seebeck  coefficient  (V  K  1 ) 

h 

convection  heat  transfer  coefficient  (W  m  2  K  ') 

ria 

absorptivity  of  the  collector  coating 

I 

electric  current  (A) 

riopt 

optical  efficiency  of  the  Fresnel  lens 

IL 

electric  current  load  vector  (A) 

P 

electrical  resistivity  (£2  m) 

J 

electric  current  density  vector  (A  m  2) 

Peon 

electrical  contact  resistivity  (£2  m  2) 

I<TT 

thermal  stiffness  matrix 

4> 

electric  scalar  potential  (V) 

K<pt 

I(V<P 

k 

Seebeck  stiffness  matrix 
electric  stiffness  matrix 
thermal  conductivity  (W  m  1  K  ’) 

cpe  vector  of  nodal  electric  potentials  (V) 

Subscripts 

C  thermoelectric  cooler 

L 

length  of  thermoelectric  element  (mm) 

N 

element  shape  function 

corn 

heat  convection  between  the  thermoelectric  elements 

q" 

heat  flux  (W  m  2) 

and  the  environment 

Qs 

solar  irradiance  (W  m  2) 

G 

thermoelectric  generator 

Q 

heat  transfer  rate  (W) 

H 

hot  side 

Qi„ 

input  energy  of  solar  thermoelectric  generator  (W) 

L 

cold  side 

0.TE 

input  energy  per  thermoelectric  couple  (W) 

N-type  thermoelectric  element 

Qt 

R 

T 

thermal  load  vector  (W) 
electrical  resistance  (£2) 
temperature  (°C) 

P 

P-type  thermoelectric  element 

the  generation  mode,  Jang  et  al.  [16]  analyzed  the  geometric  effect 
of  thermoelectric  elements  on  micro-TEG  performance  by  means  of 
a  finite  element  method.  Their  results  showed  that  there  was  an 
optimal  length  of  thermoelectric  elements  to  achieve  the  highest 
power.  Additionally,  a  higher  efficiency  could  be  obtained  with  a 
greater  length  of  thermoelectric  element.  For  the  cooling  mode, 
Lee  and  Kim  [17]  used  a  numerical  method  to  investigate  the  cool¬ 
ing  performance  of  a  micro-TEC.  They  reported  that  the  cooling 
rate  increased  and  the  maximum  COP  decreased  when  the  thick¬ 
ness  of  the  thermoelectric  element  decreased. 

As  far  as  the  practical  thermoelectric  system  is  concerned,  in  fact, 
there  exists  an  undesired  effect  of  electrical  contact  resistance  con¬ 
tributed  from  both  interfaces  and  interconnects  [18],  Min  and  Rowe 
[19]  reported  that  the  effect  of  contact  resistance  on  the  COP  of  a  TEC 
became  significant  when  the  length  of  thermoelectric  element  was 
relatively  short.  Another  unfavorable  effect  regarding  geometry  is 
heat  loss  to  the  environment  from  TEG  surfaces.  The  effect  of  heat 
loss  from  thermoelectric  elements  to  the  ambient  environment  on 
TEG  performance  has  been  explored  in  recent  research.  For  instance, 
a  TEG  system  combined  with  parallel-plate  heat  exchangers  was 
constructed  by  Niu  et  al.  [20]  where  a  hot  fluid  and  a  cold  fluid 
passed  through  the  hot  side  and  the  cold  side  of  the  TEG,  respec¬ 
tively.  Their  experimental  results  indicated  that  heat  loss  from  the 
TEG  to  the  environment  increased  markedly  when  the  inlet  fluid 
temperature  at  the  hot  side  was  lifted.  The  theoretical  analysis  of 
Xiao  etal.  [21]  revealed  that  the  existence  of  heat  loss  causedalarge 
exergy  loss  in  a  TEG  system  and  this  effect  should  be  taken  into 
account  in  the  analysis.  However,  the  geometric  design  of  the  ther¬ 
moelectric  system  was  not  considered  in  their  study. 

From  the  review  of  the  above  literature,  it  is  evident  that  spec¬ 
ifying  optimal  geometric  design  parameters  may  be  a  promising 
method  of  improving  or  maximizing  the  performance  of  thermo¬ 
electric  systems.  However,  to  the  authors’  knowledge,  very  little 
research  has  been  performed  on  the  geometric  design  of  thermo¬ 
electric  systems  under  the  effects  of  heat  loss  and  electrical  contact 
resistance,  especially  in  an  integrated  TEG-TEC  system.  In  order  to 
provide  a  useful  insight  into  the  importance  of  geometric  design 
for  improving  the  performance  of  an  integrated  TEG-TEC  system, 


a  numerical  method  is  developed  to  model  an  integrated  system 
and  predict  the  performance  of  the  system.  The  effects  of  heat  loss 
and  electrical  contact  resistance  on  the  performance  are  taken  into 
account.  The  effects  of  two  different  boundary  conditions  due  to 
altered  operating  conditions  are  also  considered. 


2.  Methodology 

2.J.  Physical  model  and  assumptions 

A  schematic  of  the  integrated  system  for  study  is  shown  in 
Fig.  la  where  a  TEC  and  a  TEG  are  included  in  the  system.  In  the 
system,  the  TEG  absorbs  heat  QhiC  from  the  heat  source  at  the 
hot  side  and  liberates  heat  Q,iC,  to  the  heat  sink  at  the  cold  side, 
thereby  generating  electric  current  through  the  Seebeck  effect 
[22],  Then,  the  current  is  directly  used  to  power  the  TEC,  which 
absorbs  heat  Q/  C  from  the  refrigerated  object  at  the  cold  side  and 
dissipates  heat  Qh.c  to  the  heat  sink  at  the  hot  side  based  on  the 
Peltier  effect  [23],  In  order  to  investigate  the  influence  of  heat  loss 
on  the  performance  of  the  integrated  system,  a  heat  convection 
process  between  the  thermoelectric  elements  and  its  environment 
is  taken  into  account.  The  radiation  heat  transfer  is  relatively  small 
at  low-temperature  conditions  [21],  so  it  is  neglected  in  this  study. 
In  Fig.  la,  Qconv.c  and  Qc0nv,c  denote  the  convective  heat  flow  rates 
along  the  lateral  surfaces  of  the  TEG  and  the  TEC,  respectively. 
For  the  sake  of  simplicity,  the  following  assumptions  are  adopted. 

(a)  The  integrated  system  is  in  steady-state. 

(b)  An  identical  model  is  used  for  the  TEC  and  the  TEG;  there¬ 
fore,  they  have  the  same  configurations  and  material  proper¬ 
ties,  except  for  the  Seebeck  coefficients  which  are  positive 
and  negative  in  the  p-type  and  n-type  elements, 
respectively. 

(c)  Material  properties  of  the  thermoelectric  elements  are  tem¬ 
perature-dependent. 

(d)  The  thermoelectric  elements  are  connected  electrically  in 
series  and  thermally  in  parallel. 


(a) 
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|  Qh,c 


(b) 


Fig.  1.  Schematics  of  (a)  the  integrated  TEG-TEC  system  for  investigation  and  (b)  thermoelectric  element  geometry. 


(e)  The  contact  resistivity  is  assumed  to  be  constant. 

(f)  In  the  TEG  or  the  TEC,  only  one  thermoelectric  pair,  namely, 
a  p-type  element  and  an  n-type  element,  are  considered. 

Accordingly,  the  computational  geometries  of  the  TEG  and  TEC 
are  shown  in  Fig.  lb.  Each  thermoelectric  element  is  specified  by 
its  depth  D,  width  W,  and  length  I.  The  lateral  surfaces  considering 
convective  heat  transfer  exclude  the  hot-side  and  cold-side 
surfaces. 

2.2.  Governing  equations 

For  a  steady-state  thermoelectric  model,  the  governing  equa¬ 
tions  include  the  thermal  and  electrical  fields  to  describe  the  ther¬ 
moelectric  effects  and  they  are  written  as  [24,25] 

V-(arJ)-V-(kVT)=J.£  (1) 


where  T  and  J  represent  the  absolute  temperature  and  electric  cur¬ 
rent  density  vector,  respectively;  a,  k,  and  p  are  the  Seebeck  coeffi¬ 
cient,  thermal  conductivity,  and  electrical  resistivity  of  a 
thermoelectric  element,  respectively.  E  is  the  electric  field  and  it 
can  be  derived  from  electric  scalar  potential  (j>  (E  =  -V<£).  More 
detailed  derivations  of  the  governing  equations  can  be  found  in  a 
previous  study  [25]  in  which  a  TEG  system  has  been  successfully 
modeled  by  the  numerical  method. 

2.3.  Boundary  conditions 

The  boundary  conditions  are  made  up  of  three  different  sur¬ 
faces,  including  hot-side  surfaces,  cold-side  surfaces,  and  lateral 
surfaces.  The  hot-side  and  cold-side  surfaces  of  the  TEC  and  TEG 
elements  are  treated  by  Dirichlet  conditions;  that  is,  these  surfaces’ 


temperatures  are  fixed.  Another  case  with  Neumann  conditions  for 
the  hot-side  surfaces  of  the  TEG  elements  are  considered  for  com¬ 
parison.  In  this  case,  a  heat  transfer  rate  is  given  at  the  hot-side 
surfaces  of  the  TEG  elements.  The  lateral  surfaces  of  all  the  thermo¬ 
electric  elements  are  exposed  to  the  environment  along  with  con¬ 
vective  heat  transfer  as  mentioned  earlier.  Therefore,  a  uniform 
convective  heat  transfer  coefficient  hconv  is  invoked  at  the  lateral 
surfaces.  Detailed  boundary  conditions  are  summarized  in  Table  1 . 

2.4.  Numerical  method  and  system  performance 

The  governing  equations  are  solved  by  means  of  the  commercial 
software  ANSYS  V12.0.1  in  which  a  finite  element  scheme  based  on 
the  Galerkin  method  is  used  to  discretize  the  governing  equations 
[26],  Specifically,  the  physical  scales  of  temperature  T  and  electric 
scalar  potential  <j>  over  a  computational  domain  are  approached  by 
the  finite-element  method  by  the  following  equations  [27] 

T=[N]{Te}  (3) 

4>=[N]{<pJ  (4) 

where  Te,  tpe,  and  N  are  the  vector  of  nodal  temperature,  the  vector 
of  nodal  electrical  potentials,  and  the  element  shape  function, 
respectively.  By  integrating  Eqs.  (1)  and  (2)  based  on  the  Galerkin 


A  list  of  adopted  boundary  conditions. 


Surface 

TEC 

TEG 

Hot  side 

Dirichlet  condition: 

T=Th.c 

Casel :  Dirichlet  condition: 

T=T„,C 

Case2:  Neumann  condition: 

/'  kVTdA  =  q"A  =  QHG 

Cold  side 

Dirichlet  condition:  T=  T^c 

Dirichlet  condition:  T=  \c 

Convection  condition: 

Convection  condition: 

side 

-kvr=h„nv(r-r00) 

-kVT=hconv(T-TJ 
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method,  the  differential  equations  can  be  rewritten  as  algebraic 
equations  [27] 

[s  ;.]{;•>{?} 

where  K77,  Krprp,  and  KtpT  are  the  thermal  stiffness  matrix,  electric 
stiffness  matrix,  and  Seebeck  stiffness  matrix,  respectively,  and  they  ~ 
are  defined  as  K77  =  f„VN ■  [fc]  ■  VNdV,  K w  =f„VN-  [±1  ■  VNdV.  % 

and  K,pJ  =  fvVN  ■  [i]  •  [a]  •  VNdV.  Q'  and  f  in  Eq.  ( 5 )  are  the  thermal  X 

load  vector  and  the  electric  current  load  vector,  respectively.  The  O 
three  matrices  and  two  load  vectors  are  obtained  by  the  numerical 
integration  over  the  element  volume.  Four  thermoelectric  elements 
are  considered  in  the  numerical  model,  as  shown  in  Fig.  lb.  The 
orthogonal  grid  system  is  used  for  the  thermoelectric  elements. 

The  test  of  grid  system  reveals  that  each  element  consisting  of 
5850  numerical  cells  satisfies  the  requirement  of  grid  indepen¬ 
dence;  this  number  of  cells  is  thus  adopted  for  simulations.  In  the 
calculation,  the  iterations  are  terminated  as  the  residuals  of  all 
the  equations  reach  the  convergence  criteria  of  relative  errors  less 
than  10  6.  After  the  numerical  iterations  are  implemented,  the  tem¬ 
peratures  (Te)  and  electric  potentials  (<j>e)  at  each  node  of  the 
numerical  cells  are  obtained. 

In  the  integrated  system,  two  indices  are  employed  to  evaluate 
the  performance  of  the  system.  One  is  the  cooling  power  Qj  C  which 
represents  the  heat  pumped  capacity  of  the  TEC  from  the  refriger¬ 
ated  object.  The  other  is  the  COP  of  the  system  and  it  is  expressed 
as  [11,14] 


COP  means  the  percentage  of  input  energy  (i.e.  Qh.g)  gained  for 
cooling  power  by  the  system.  The  thermoelectric  module  TEC1- 
12708  presented  by  Maneewan  and  Chindaruksa  [28]  was  chosen 
as  the  geometric  design  of  the  TEC  and  the  TEG  in  this  study,  as 
given  in  Table  2.  The  temperature-dependent  material  properties 
of  the  thermoelectric  elements,  such  as  thermal  conductivity,  See¬ 
beck  coefficient,  and  electrical  resistivity,  reported  by  Meng  et  al. 
[29]  were  adopted  and  listed  in  Table  2.  The  material  is  a  commer¬ 
cially  available  material  produced  by  Melcor;  however,  detailed 
components  in  the  material  were  not  illustrated,  perhaps  due  to 
the  commercial  know-how  involved. 

In  reviewing  past  studies,  very  little  experimental  literature  and 
available  data  have  been  published  on  the  integrated  system, 
Therefore,  theoretical  solutions  were  employed  to  validate  numer¬ 
ical  predictions  in  some  studies  [30-33].  Similarly,  the  present 
model  is  validated  by  comparing  with  the  analysis  of  Khattab 
and  Shenawy  [12],  Their  method  was  conducted  based  on  one¬ 
dimensional  (1-D)  geometry  and  the  assumption  of  ignoring  heat 
convection  between  the  thermoelectric  elements  and  their  envi¬ 
ronment;  the  formulas  are  given  by: 

Q.h,g  =  «JTl,c  +  0.5 I2Rc  +  Kc(Th, c  ~  Tyj)  (7) 


75  80  85  90 

Th,g(°C) 


Fig.  2.  A  comparison  of  one-dimensional  system  performance  between  numerical 
simulation  and  theoretical  prediction. 

Qi,c  =  Ocffqr  -  0.5 12RC  -  KC{TH,C  -  Tqc)  (8) 

,  _  <Xg(Th,g  -  TLfi)  -  0Cc(Th,c  -  Tl,c ) 

1  ~  Rg  +  Rc  a 

where  y.c  and  ac  are  the  Seebeck  coefficient  of  the  TEG  and  TEC, 
respectively;  Rc  and  Rc  represent  the  resistance  of  the  TEG  and  of 
the  TEC,  respectively;  I<c  and  Kc  are  the  thermal  conductance  of 
the  TEG  and  of  the  TEC,  respectively.  The  assumptions  of  ocG  =  «c  = 
(otp  -  o„),  Rc  =  Rc  =  ^  and  Kc  =  Kc  =  ^  ^  are  consid¬ 
ered.  Therefore,  the  theoretical  COP  can  be  obtained  from 
COP  =  QL'd Qh.g-  By  virtue  of  the  limitation  of  1-D  approach,  the 
constant  material  properties  at  T  =  300  I<  were  used  and  the  contact 
resistivity  pcon  was  set  to  be  1  x  10  10  £1  m  2  [34],  The  numerical 
and  theoretical  results  along  with  T,  Cl  =  25  °C,  THC  =  25  °C,  and 
TLc=  15  °C  are  shown  in  Fig.  2.  The  numerical  predictions  are  in 
excellent  agreement  with  the  theoretical  results,  so  the  developed 
model  is  reliable  for  analyzing  the  performance  of  the  integrated 
system. 

2.5.  Investigation  basis 

The  element  lengths  of  TEG  and  TEC  ranging  from  0.6  to  2.4  mm 
are  considered  to  account  for  the  influence  of  geometry  on  system 
performance.  Two  different  boundary  conditions  at  the  hot  surface 
of  TEG,  with  one  the  Neumann  condition  (i.e.  a  temperature  is 
given)  and  the  other  with  the  Dirichlet  condition  (i.e.  a  heat  trans¬ 
fer  rate  is  given),  are  considered.  The  geometric  values  of  base  case 
and  the  operating  (boundary)  conditions  are  shown  in  Tables  2  and 


material  properties  of  thermoelectric  elements. 


Geometry  ofTECl-12708  [28] 

D  =  1.4  mm 
1  =  1.2  mm 


Material  properties  [29] 

k  =  (62605  -  277.7T+  0.4131T2)  x  lO^Wm-1  K  1 
p  =  (5112  +  163.4T+  0.6279T2)  x  lO”10  fim 
a  =  (22224  +  930.6T  -  0.9905T2)  x  109  V  K  1 


Table  3 

Operating  (boundary  condition)  conditions  of  study. 


Surface  TEC 

Hot  side  Dirichlet  condition: 

T„,c=  25  °C 


Cold  side  Dirichlet  condition: 
Tl,c=  15  °C 

Lateral  Convection  condition: 

side  h„„v  =  20Wm-2K-1 

Tx,  =  25  °C 


TEG 


Casel:  Dirichlet  condition: 

Th.c=  100  X 

Case2:  Neumann  condition: 

Qh.g  =  0.5  W 

Dirichlet  condition:  Tuc  =  25  °C  °C 

Convection  condition 
h„  =  20Wm-2r' 

T„.  =  25  °C 


570 


W.-H. 


:  al./ Energy  Conversion  and  Management  87  (2014)  566-575 


(b)  0105 

I  ;  :  :  ;  . 

0.095 

O  0.09 
V 

0.085 

7  v  "  * 

-  Ideal  case  # 

0  08 

-  -  Ideal  case 

_  ^  Pc°n-2xl010  Qm2  T 

- T - PcZ=4xlO'10  (1m2 

0  075 

-  - T - p>0-10  Qm1 

0  0  10  20  30  40  50  60 

hc„„v  (W/m2K)  hconv(W/m2K) 


Fig.  3.  Distributions  of  (a)  hot  surface  temperature  of  the  TEG  and  (b)  electric  FiS-  4‘  Distributions  of  (a)  cooling  power  and  (b)  COP  (base  case), 

current  of  the  system  (base  case). 


3,  respectively.  The  solar  thermoelectric  generator,  which  uses 
concentrated  solar  radiation  as  a  heat  source,  has  received  much 
attention  lately  [35],  The  input  energy  (Qjn)  of  a  thermal-concen¬ 
trated  solar  thermoelectric  generator  can  be  represented  by  [36,37] 


aforementioned  boundary  condition.  The  contact  resistivity  and 
convection  heat  transfer  coefficient  of  the  base  case  are  set  as 
1  x  10  lofim-2  [34]  and  20  Wm  2K  1  [21,38],  respectively. 

3.  Results  and  discussion 


Qin  =  qs  x  Cg  xAc  x  rjopt  x  r,a  (10) 

where  qs,  Cg,  Ac,  r]opt,  and  qa  are  the  solar  irradiance,  the  concentra¬ 
tion  ratio,  the  cross-sectional  area  of  the  collector,  the  optical  effi¬ 
ciency  of  the  Fresnel  lens,  and  the  absorptivity  of  the  collector 
coating,  respectively.  The  cross-sectional  area  of  the  collector  is 
assumed  to  be  equivalent  to  the  cross-sectional  area  of  the  thermo¬ 
electric  module.  With  the  conditions  of  qs  =  900  W  m  2,  Cg  =  50, 
qopt  =  85,  and  qa  =  0.9  [36]  as  well  as  Ac  =  40  x  40  mm2  for  the  ther¬ 
moelectric  module  TEC1-12708  [28],  the  input  energy  is 

Qin  =  900  x  50  x  (40  x  40  x  10“6)  x  0.85  x  0.9  =  55.08  W  (11) 
TEC1-12708  consists  of  127  thermoelectric  couples  [28],  so  the 
input  energy  per  couple  (Qtb)  is  Qje  =  55.08/127  a;  0.43  W  and  this 
value  is  close  to  0.5  W.  For  this  reason,  the  boundary  condition  of 
Qh,c  =  0.5  W  is  selected  in  this  work.  In  addition,  the  maximum 
temperature  of  the  TEG  does  not  exceed  the  operating  temperature 
limit  of  the  thermoelectric  module  in  all  simulations  with  the 


3.1.  Influences  of  heat  convection  and  contact  resistance  on  system 
performance 

The  influences  of  convection  heat  transfer  coefficient  and  contact 
resistivity  on  system  performance  at  the  condition  of  Qh.g  =  0.5  W 
are  shown  in  Figs.  3  and  4.  The  convection  heat  transfer  coefficient 
and  contact  resistivity  are  in  the  ranges  of  0-50  W  m  2  K  1  [21,37] 
andl  x  10  10-4  x  10  10  Dm  2  [18,34],  respectively.  The  ideal  case 
disregarding  heat  convection  and  contact  resistance  is  also  included 
for  comparison.  It  can  be  seen  that  the  hot  surface  temperature  of  the 
TEG  rises  with  increasing  pcon  but  decreases  with  hconv  (Fig.  3a).  The 
energy  balance  at  the  TEG  hot  surface  is  expressed  as  [7,21  ] 

Qh,G  =  Qpe  +  Qf  -  Qj  +  Qconvfi  (12) 

where  Qpe,  Qr,  0/,  and  Qconv.c  designate  the  heat  transfer  rate  due  to 
the  Peltier  effect,  Fourier  conduction,  Joule  heating,  and  heat  con¬ 
vection,  respectively.  The  contact  resistance  can  be  conceived  as 
an  additional  internal  resistance  of  thermoelectric  module  [39]; 
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hence,  an  increase  in  pcon  raises  the  Joule  heat.  In  order  to  keep 
energy  balance,  the  increased  Joule  heat  may  lead  to  the  increase 
of  Peltier  heat  (i.e.,  Eq.  (12)).  Accordingly,  the  hot  surface  tempera¬ 
ture  of  the  TEG  increases  because  the  Peltier  heat  at  the  hot  surface 
is  expressed  as  Qpe  =  (ap  -  an)/TH,G  [7,21],  Similarly,  increasing  con¬ 
vection  heat  transfer  coefficient  abates  the  Peltier  heat  due  to  the 
fixed  Qh.g,  thereby  decreasing  the  hot  surface  temperature  of  the 
TEG.  The  electric  current  of  the  integrated  system  diminishes  as  pcon 
goes  up  (Fig.  3b).  This  is  attributed  to  the  increased  total  resistance 
of  the  integrated  system.  The  electric  current  also  declines  with 
increasing  hconv,  and  this  reduces  the  cooling  power  (Fig.  4a).  This 
is  because  the  heat  pumped  at  the  TEC  cold  side  through  the  Peltier 
effect  (expressed  as  ( ap  -  an)/T/  c  [30])  is  the  dominant  mechanism 
[24],  A  similar  trend  in  COP  under  the  condition  of  fixed  Qh.g  is  also 
exhibited  (Fig.  4b).  Compared  to  the  ideal  case,  the  maximum 
reduction  in  cooling  power  (or  COP)  is  18.11%  which  occurs  at 
pcon  =  4  x  10  10  Q  m  2  and  hconv  =  50  W  m  2  K 

3.2.  Influence  of  TEC  element  length  at  a  given  hot  surface  temperature 
ofTEC 

Two  kinds  of  boundary  condition  at  the  hot  surface  of  TEG  are 
employed  to  analyze  the  influence  of  TEG  element  length  on  sys¬ 
tem  performance.  Meanwhile,  three  different  models  are  regarded 


for  comparison:  Model  1  (without  contact  resistance  and  heat  con¬ 
vection),  Model  2  (with  contact  resistance  but  without  heat  con¬ 
vection),  and  Model  3  (with  both  contact  resistance  and  heat 
convection).  With  the  condition  of  TH,C  =  100  °C,  the  electric  cur¬ 
rent  of  the  integrated  system  decreases  as  the  TEG  length  increases 
(Fig.  5a),  resulting  from  the  increased  internal  resistance  of  the 
TEG.  This  phenomenon  can  also  be  found  from  the  theoretical  for¬ 
mula  of  Khattab  and  Shenawy  [12]  (i.e.,  Eq.  (9)).  The  decrease  of 
electric  current  also  abates  the  Peltier  heat  at  the  TEG  hot  surface. 
Therefore,  Qh.g  is  lessened  as  the  TEG  element  length  increases 
(Fig.  5b).  In  Model  2,  because  the  contact  resistance  increases  the 
total  resistance  of  the  integrated  system,  the  electric  current  of 
the  system  decreases  (Fig.  5a).  Furthermore,  the  decreased  electric 
current  and  the  increased  resistance  make  the  Joule  heat  become 
negligible.  This  is  the  reason  why  the  heat  input  values  of  Model 
2  are  close  to  those  of  Model  1  (Fig.  5b).  Seeing  that  the  two  surface 
temperatures  of  the  TEG  are  fixed,  the  additional  heat  convection 
considered  in  Model  3  plays  no  part  in  changing  the  electric  cur¬ 
rent  when  compared  with  Model  2  (Fig.  5a);  however,  a  bit  more 
heat  input  of  Qh.g  is  required  in  Model  3. 

The  Peltier  heat  at  the  cold  side  of  the  TEC  is  the  dominant 
mechanism  of  cooling  power  (Q^c)  and  it  is  proportional  to  the 
electric  current;  therefore,  Qic  decreases  with  increasing  TEG 
length  (Fig.  6a).  When  the  effect  of  contact  resistance  on  the 
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reduction  percentage  of  Qr,c  is  examined,  a  minimum  value  devel¬ 
ops  at  Lg  =  1 .8  mm.  If  both  the  heat  convection  and  contact  resis¬ 
tance  are  simultaneously  considered,  the  maximum  reduction 
percentage  of  Ql,c  takes  place  at  Lc  =  2.4  mm.  As  a  whole,  both 
Qh.c  and  Qf.,c  are  the  decreasing  functions  of  TEG  element  length 
(Figs.  5b  and  6a).  It  is  worthy  of  note  that  there  exists  an  optimum 
Lc  corresponding  to  the  maximum  COP  (Fig.  6b).  For  example,  the 
optimum  TEG  length  in  Model  3  is  0.9  mm.  With  regard  to  the 
reduction  percentage  of  COP,  Fig.  6b  depicts  that  the  effect  of  heat 
convection  on  COP  tends  to  become  more  pronounced  than  that  of 
contact  resistance  when  the  element  length  goes  up.  Within  the 
investigated  range  of  TEG  length,  the  reduction  percentage  of 
COP  is  between  3.47%  and  10.53%. 


3.3.  Influence  of  TEG  element  length  at  a  given  heat  transfer  rate  of 
TEG 

Figs.  7  and  8  demonstrate  the  system  performance  under  a 
given  heat  transfer  rate  of  QhiC  =  0.5  W  at  the  hot  surface  of  the 
TEG,  whereas  its  cold  surface  temperature  is  fixed.  Contrary  to 
Fig.  5  in  which  the  hot  surface  temperature  is  fixed,  an  increase 
in  element  length  raises  the  TEG  hot  surface  temperature 


(Fig.  7a),  stemming  from  the  increased  thermal  resistance  of  the 
TEG  [40],  In  view  of  the  enlarged  temperature  difference,  the  elec¬ 
tric  current  of  the  system  grows  with  increasing  TEG  element 
length.  Eq.  (9)  also  elucidates  the  preceding  behavior.  The  electric 
current  of  the  system  increases  slowly  when  the  length  is  large  to  a 
certain  extent  (Fig.  7b).  The  output  power  of  the  TEG  is  0.018  W 
when  its  hot-side  temperature  reaches  about  145  °C.  This  power 
is  very  close  to  the  experimental  result  of  Niu  et  al.  [20]  where 
the  TEG  power  was  0.0183  W  at  the  same  hot  surface  temperature. 
As  illustrated  earlier,  the  electric  current  directly  affects  the  Peltier 
heat  at  the  cold  surface  of  the  TEC;  hence  increasing  TEG  element 
length  increases  the  cooling  power  and  COP  (Fig.  8)  and  their  pro¬ 
files  are  similar  to  those  of  electric  current.  It  should  be  pointed  out 
that  the  case  of  Lc  =  2.4  mm  along  with  Qh.c  =  0.5  W  gives  the  best 
performance  of  the  system  (Fig.  8),  but  the  lowest  performance  is 
exhibited  at  the  condition  of  THC  =  100  °C  (Fig.  6).  With  attention 
paid  to  the  reduction  percentages  of  cooling  power  and  COP, 
Fig.  8  reveals  that  the  reduction  due  to  the  contact  resistance 
declines  when  the  element  length  increases,  and  the  impact  of  con¬ 
vection  on  the  reduction  of  cooling  power  or  COP  is  more  than  the 
contact  resistance.  This  is  the  reason  that  the  difference  between 
the  curves  of  Models  2  and  3  is  more  notable  than  that  between 
Models  1  and  2.  Overall,  the  reduction  percentage  of  COP  caused 
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by  the  heat  convection  and  contact  resistance  is  in  the  range  of 
6.22-12.45%. 

3.4.  Influence  of  TEC  element  length  on  system  performance 

The  influence  of  TEC  element  length  on  system  performance  at 
the  condition  of  THjC  =  100  °C  is  examined  in  Figs.  9  and  10  where 
the  three  different  models  are  also  included.  The  internal  resis¬ 
tance  of  TEC  increases  when  its  element  length  rises,  thereby 
diminishing  the  electric  current  of  the  system  (Fig.  9a).  This  further 
lessens  the  Peltier  heat  at  the  TEG  hot  surface  (i.e.  Qpe  =  (ocp  -  an)~ 
TTh.c)  so  as  to  decrease  the  heat  input  into  the  TEG  (Fig.  9b).  It  is 
worth  noting  that  an  optimum  Ic  corresponding  to  the  maximum 
cooling  power  and  COP  can  be  identified  (Fig.  10).  For  a  shorter 
TEC,  the  importance  of  Peltier  heat  at  the  TEC  cold  surface  precedes 
the  Joule  heat,  as  consequences  of  larger  electric  current  and  smal¬ 
ler  internal  resistance.  However,  heat  diffusing  back  to  the  cold 
surface  of  the  TEC  by  the  Fourier  conduction  is  pronounced 
because  the  total  conductance  of  a  thermoelectric  couple  is  inver¬ 
sely  proportional  to  the  element  length  [14].  A  lower  cooling 
power  is  thus  featured  by  a  shorter  TEC.  When  the  TEC  becomes 
taller,  the  Fourier  conduction  is  weakened  so  the  cooling  power 
increases.  When  Lc  is  large  to  a  certain  extent,  the  Peltier  heat 


decreases  because  of  smaller  electric  current.  Once  the  Peltier  heat 
is  down  to  a  level  insufficient  to  counteract  the  Joule  heat  and  Fou¬ 
rier  conduction  heat,  the  cooling  power  decreases.  Unlike  the 
behavior  shown  in  Fig.  5b,  there  exist  perceptible  differences  in 
Qh.c  among  the  three  curves  shown  in  Fig.  9b,  implying  that  the 
existences  of  contact  resistance  and  convection  affect  the  heat 
input  to  a  small  extent,  especially  for  the  latter  where  more  heat 
inputted  to  the  TEG  is  observed.  In  Fig.  10,  the  largest  total  reduc¬ 
tion  percentages  of  cooling  power  and  COP  are  7.10%  and  7.65%, 
respectively,  where  the  TEC  length  (Lc)  is  0.6  mm. 

Upon  inspection  of  system  performance  at  CW  =  0.5  W  under 
the  situation  of  varied  TEC  length,  the  electric  current  of  the  system 
decreases  with  increasing  TEC  length  (Fig.  11a),  resulting  from  the 
increased  internal  resistance  of  the  TEC.  Furthermore,  the  reduced 
current  affects  the  energy  balance  at  the  TEG  hot  surface.  On 
account  of  fixed  CW  and  in  order  to  keep  the  energy  balance  at 
the  TEG  hot  surface,  the  hot  surface  temperature  of  the  TEG 
increases  (Fig.  lib).  However,  the  temperature  variation  versus 
the  length  is  not  as  significant  as  that  shown  in  Fig.  7a.  It  is  also 
noted  that  the  contact  resistance  lifts  the  TEG  hot  surface  temper¬ 
ature.  The  characteristics  of  the  maximum  cooling  power  and  COP 
are  observed  in  Fig.  12  where  the  optimum  TEC  length  is  1.8  mm. 
The  foregoing  results  are  similar  to  those  with  fixed  THC  (Fig.  10). 
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Fig.  11.  Distributions  of  (a)  electric  current  of  the  system  and  (b)  hot  surface 
temperature  of  the  TEG  (Qh,g  =  0.5  W). 


Fig.  12.  Distributions  of  (a)  cooling  power  and  (b)  COP  and  their  reduction 
percentages  (Qh.c  =  0.5  W). 


This  can  be  explained  by  the  same  trend  in  electric  current,  even 
though  the  boundary  conditions  at  the  TEG  hot  surface  are  differ¬ 
ent  from  each  other.  The  negative  effects  of  contact  resistance 
and  heat  convection  on  system  performance  can  be  found  in 
Figs.  11  and  12.  The  largest  total  reduction  percentages  of  cooling 
power  and  COP  are  obtained  at  Lc  =  0.6  mm  where  their  values  are 
the  same  as  18.67%. 

The  characteristics  of  cooling  power  and  COP  in  Model  3  in 
accordance  with  the  selection  of  boundary  condition  and  the  vari¬ 
ation  of  TEG  or  TEC  length  are  summarized  in  Table  4.  Apparently, 
the  characteristics  of  cooling  power  and  COP  under  the  condition 
of  varying  TEG  length  are  fairly  different  from  those  of  varying 
TEC  length.  The  boundary  conditions  play  an  important  role  in 
determining  the  system  performance  when  the  TEG  length  is 
altered,  whereas  an  optimum  TEC  length  can  be  obtained,  regard¬ 
less  of  the  Dirichlet  or  Neumann  condition  adopted.  Overall,  the 
predictions  reveal  that  the  cooling  power  of  the  TEC  is  in  the  range 
of  approximately  10-100  mW,  which  is  by  far  smaller  than  the 
requirement  of  practical  application.  Accordingly,  a  feasible  way 
for  the  application  of  the  integrated  system  is  that  a  TEC  can  couple 
with  a  number  of  TEGs  in  series  to  intensify  the  cooling  power  of 
the  TEC.  The  aforementioned  system  deserves  further  investigation 
in  the  future. 


Summary  of  system  performance  of  Model  3  at  various  operating  conditions. 

Cooling  power  COP 

Increasing  TEC  length 

Dirichlet  condition3  Decrease  Optimum  at  Ic  =  0.9  mm 

Neumann  condition13  Increase  Increase 

Increasing  TEC  length 

Dirichlet  condition3  Optimum  at  Lc=  1.2  mm  Optimum  at  Lc=  1.2  mm 
Neumann  condition13  Optimum  at  Lc=  1.8  mm  Optimum  at  Lc=  1.8  mm 


3  TH,c  =  100  °C. 
b  Qh,g  =  0.5W. 


4.  Conclusions 

The  physical  phenomena  of  an  integrated  TEG-TEC  system  have 
been  analyzed  through  a  numerical  method.  Particular  attention  is 
paid  to  the  influences  of  geometric  sizes  of  TEG  and  TEC  and 
boundary  conditions  upon  the  performance  of  the  integrated  sys¬ 
tem.  The  effects  of  contact  resistance  and  heat  convection  on  the 
performance  have  also  been  examined.  The  predictions  suggest 
that  the  cooling  power  and  coefficient  of  performance  (COP)  of 
the  system  are  reduced  when  the  contact  resistance  and  heat  con¬ 
vection  are  taken  into  account.  For  a  given  hot  surface  temperature 
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of  TEG,  the  higher  the  TEG  length,  the  larger  the  reduction  percent¬ 
ages  of  cooling  power  and  COP.  For  the  other  cases,  an  optimum 
length  for  minimizing  the  reduction  percentage  can  be  found. 
However,  the  minimum  reduction  percentage  does  not  always 
respond  to  the  maximum  system  performance.  When  the  lengths 
of  TEG  and  TEC  are  individually  altered,  the  system  performances 
are  reduced  up  to  12.45%  and  18.67%,  respectively.  Considering 
the  impacts  of  geometric  size  and  boundary  conditions  on  the  per¬ 
formance,  the  cooling  power  and  COP  of  the  system  are  signifi¬ 
cantly  affected  by  the  boundary  conditions  when  the  TEG  length 
is  changed.  In  contrast,  an  optimum  length  corresponding  to  the 
maximum  cooling  power  and  COP  is  exhibited  if  the  TEC  length 
is  varied,  irrespective  of  whether  the  boundary  condition  of  fixed 
temperature  (Dirichlet)  or  heat  transfer  rate  (Neumann)  is  given. 
The  predictions  suggest  that  the  cooling  power  of  the  TEC  in  the 
integrated  system  is  in  the  range  of  approximately  10-100mW. 
These  values  are  fairly  small  for  practical  applications.  A  feasible 
way  for  the  application  of  the  integrated  system  is  the  coupling 
of  a  TEC  with  a  number  of  TEGs  in  series  to  intensify  the  cooling 
power  of  the  TEC.  The  present  study  has  provided  a  useful  insight 
into  the  design  of  integrated  TEG-TEC  systems  which  enables  us  to 
refrigerate  an  object  by  using  waste  heat  as  a  power  source. 
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